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Introduction. Two prominent methods to treat non-perturbative Yang-Mills theory will be addressed in this talk: 
the numerical simulation of lattice gauge theory (LGT) and the approach by means of the Dyson- Schwinger equa- 
tions (DSE). While LGT covers all non-perturbative effects and, in particular, bears witness of quark confinement 
(see e.g. fllf) , simulations including dynamical quarks are cumbersome despite the recent successes by improved al- 
gorithms B and the increase of computational power. At the present stage, systems at finite baryon densities are 
hardly accessible in the realistic case of an SU(3) gauge group |J (for recent successes see Q). By contrast, the 
DSE approach can be easily extended for an investigation of quark physics even at finite baryon densities 0. 
Unfortunately, the DSE approach requires a truncation of the infinite tower of equations, and this approximation is 
difficult to improve systematically. In addition, the DSE approach needs gauge fixing which is obscured by Gribov 
copies. Whether the standard Faddeev-Popov method of gauge fixing is appropriate in non-perturbative studies, is 
still under debate ||]. 

In order to merge the advantages of both approaches to low energy Yang-Mills theory, I will firstly address the gluon 
propagator of pure SU(2) lattice gauge theory in Landau gauge. The result can then be compared with the one 
provided by the solution of the coupled ghost-gluon Dyson equation @JOjl3]. This will allow us to estimate the 
soundness of the truncations introduced to solve the equations (e.g. vertex ansatz, angular approximation). Secondly, 
the gluon propagator is an one essential ingredient of the quark DSE. Two options are obvious: a parameterization of 
the lattice result for the gluon propagator enters the quark DSE. The corresponding solution of this equation provides 
informations on hadronic observables in quenched approximation i.e. the backreaction of the quarks on the gluonic 
Greenfunctions is neglected. Once the reliability of the DSE approach to the ghost gluon system has been tested, 
the second option is to solve a truncated set of coupled ghost-gluon-quark DSEs, thereby, challenging the quenched 
approximation. 

In my talk, I will focus on the gluon propagator as inferred from the lattice calculation, and I will concentrate on 
the information on quark confinement which might be encoded in the gluon propagator. High accuracy data for the 
latter are obtained by a numerical method superior to existing techniques. Further informations and details of the 
numerical method will be presented in a forthcoming publication. 

Lattice definition of the gluon field. Before identifying the gluonic degrees of freedom in the lattice formulation, 
I briefly recall the definition of the gluon field in continuum Yang-Mills theory. 

The starting point for constructing Yang- Mills theories is the transformation property of the matter fields. In the case 
of an SU(2) gauge theory, we demand invariance under local SU(2) (say color) transformations of the quark fields 

q'{x) = Q{x) q{x) , Q(x) £ SU(2) . (1) 

In order to construct a gauge invariant kinetic term for the quark fields, one defines the gauge covariant derivative 
Dp := 0^ + iA^t a , where t a are the generators of the SU(2) gauge group. Per definitionem, this covariant derivative 
homogeneously transforms under gauge transformations, 

D'^q'ix) = Q(x)Dfa)q(x) (2) 

if the gluon fields transforms as 

Affix) = O ab (x) A b fa) + r S O ec (x) dfa^ , (3) 
O ab (x) := 2tr{n(x)t a n' ! (x)t b } , O ab (x) e SO(3) . (4) 

The crucial observation is that the gluon fields transform according to the adjoint representation while the matter 
fields are defined in the fundamental representation. 

Let us compare these definitions of fields with the ones in LGT. In LGT, a discretization of space-time with a lattice 
spacing a is instrumental. The 'actors' of the theory are SU(2) matrices C/ M (x) which are associated with the links of 
the lattice. These links transform under gauge transformations as 

Ufa) = n(x)ufa)n^(x + n) n(x)€SU(2). (5) 

For comparison with the ab initio continuum formulation, I also introduce the adjoint links 

Uf(x):=2tr{Ufa)t a Ul(x)t b }, (6) 
Ufa) := 0(x) Ufa) O t (x) , 0(x) € SO(3) , (7) 
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where 0(x) was defined in (||). 

In order to define the gluonic fields from lattice configurations, I exploit the behavior of the (continuum) gluon fields 
under gauge transformations (see (^)), and identify the lattice gluon fields A!t{x) as algebra fields of the adjoint 
representation, i.e. 



exp{e f A f (x) a} 



-i ab 



( e /)-. = e a/C ; (8) 



where the total anti-symmetric tensor e abc acts as generator for the SU(2) adjoint representation, and where a denotes 
the lattice spacing. 

For later use, it is convenient to have an explicit formula for the (lattice) gluon fields A®(x) in terms of the SU(2) 
link variables UAx). Usually, these links are given in terms of four vectors of unit length 

Up(x) = u°(z) + f, K(x)] 2 + [u^(x)] 2 = 1, (9) 

where r b are the Pauli matrices. Using these variables, a straightforward calculation yields 

A h l (x)a + Ola 2 ) — u^(x)u b ^(x) , without summation over fj, . (10) 
I point out that ([!(]) is a novel definition of the lattice gluon field. 

Finally, I illustrate the definition ([!(]). Let us assume that we have exploited the gauge degrees of freedom (see (|])) 
to bring the SU(2) link elements U^{x) as close as possible to the unit matrix, 

fl(x) : tT %( x ) -> max • ( n ) 

In this gauge, I decompose the link variables by 

U'Jx) = Z^x) exJiA b Jx)t b a) , (12) 



where A b (x) is implicitly defined by (g) and Z^(x) 6 { — 1,+1}. Indeed, the lattice gluon fields ( [Tci|) do not change 

when U^{x) — » (-l)U^(x). Hence, the fields A b (x) are relegated to the SO(3) coset space. I here propose to 
disentangle the information carried by the center and coset parts of the link variables by studying the Z^(x) and 
A b (x) correlation functions independently. I stress, however, that in Landau gauge (|TT| ) the role of the Z fl (x) is 
de-emphasized (Z^(x) — ► 1). In particular, I do not expect a vastly different gluon propagator when other (more 



standard) definitions of the lattice gluon fields are used 113 14 



Gauge fixing. In the continuum formulation, calculations employing gauge fixed Yang-Mills theory use only gauged 
gluon fields which satisfies the gauge condition, e.g. 

d^A'fix) =0 (Landau gauge) , (13) 

and rely on the assumption that the Faddeev-Popov determinant corrects the probabilistic weight in an appropriate 
way. This assumption is true if the gauge condition picks a unique solution Q(x) of (|l^) for a given field A^(x). 
Unfortunately, the Landau gauge condition generically admits several solutions depending on the "background field" 
AAx) which is the subject of gauge fixing (Gribov problem). Thus, the above assumption seems not always be 
justified ||. Further restrictions on the space of possible solutions are required [HJ. 

Let us contrast the continuum gauge fixing with its lattice analog. In a first step, link configurations U^(x) are 
generated by means of the gauge invariant action without any bias to a gauge condition. In a second step, the gauge- 
fixed ensemble is obtained by adjusting the gauge matrices Q(x) (see (j^)) until the gauged link ensemble satisfies the 
gauge condition. This procedure obviously guarantees the correct probabilistic measure for the gauged configurations, 
and gauge invariant quantities which are calculated with the gauged configurations evidently coincide with the ones 
obtained from un- fixed configurations. However, the Gribov problem re-appears as the problem of finding "the name 
of the gauge". Let me illustrate this last point: The naive Landau gauge condition for the lattice gluon field, i.e. 
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d,A'\x) = (14) 



is satisfied if we seek an extremum of the variational condition (|llj). If we restrict the variety of solutions fl(x) which 
extremize (fll) to those solutions which maximize the functional (vA\), we confine ourselves to the case where the 



Faddeev-Popov matrix is positive semi-definite. The fraction of the configuration space of gauge fixed fields A 1 ^ is 
said to lie within the first Gribov horizon. A numerical algorithm which obtains the gauge transformation matrices 
n(x) from the condition (|ll|) still samples a particular set of local maxima. Different algorithms might differ in the 
subset of chosen local maxima, and, hence, implement different gauges. A conceptual solution to the problem is to 
restrict the configuration space of gauge fixed fields A'^ to the so-called fundamental modular region. In the lattice 
simulation, this amounts to picking the global maximum of the variational condition (pd[). In practice, finding the 
global maximum is a highly non-trivial task. Here, I adopt two extreme cases of gauge fixing: firstly, I will study the 
gluon propagator of the gauge where an iteration over-relaxation algorithm almost randomly averages over the local 
maxima of the variational condition (|TT|). This result will then be compared with the gluon propagator of a gauge 
where a simulated annealing algorithm searches for the global maximum. It will turn out that the gluon propagators 
of both gauges agree within statistical error bars. 

The numerical simulation: The link configurations are generated using the Wilson action. I refrain from using a 
"perfect action" since I am interested in the gluon propagator in the full momentum range; simulations using perfect 
actions recover a good deal of continuum physics at finite values of the lattice spacing at the cost of a non-local action. 
For practical simulations, perfect actions are truncated which is poor approximation at high energies where the full 
non-locality of the action must come into play. 

Here, calculations were performed using a 16 3 x 32 lattice. The dependence of the lattice spacing on j3 (renormaliza- 
tion), i.e. 



aa 2 ([3) = 0.12 cxpj-^j-(/3- 2.3) j, a := (440 MeV) 2 , (15) 

is appropriate for [3 E [2.1,2.6] for the achieved numerical accuracy. 

Once gauge-fixed ensembles are obtained by implementing a variational gauge condition (see discussion of previous 
section), the gluon propagator is calculated using 

D*(x-y) = (Al{x)A b v {y)) MC , (16) 

where the Monte-Carlo average is taken over 200 properly thermalized gauge configurations. Of particular interest is 
the gluonic form factor F(p 2 ) which is defined by 

D(P 2 ) = J D™{x) exp{ipx} d A x , D(p 2 ) = ^ . (17) 

Since in Landau gauge the propagator is diagonal in color and transversal in Lorentz space, the form factor F(p 2 ) 
contains the full information. 



TABLE I. Simulation parameters: L = 32a: lattice extension, A: UV cutoff 



p 


2.1 


2.2 


2.3 


2.4 


2.5 


L [fm] 


8.6 


6.6 


5.0 


3.8 


2.9 


A [GeV] 


2.3 


3.0 


4.0 


5.2 


6.8 
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Results I: Landau gauge The lattice momentum in units of the lattice spacing is given by 

p x a = 2 sin ^fi^j , x = l..A (18) 

where n x is an integer which numbers the Matsubara frequency and N is the number of lattice points (in ^-direction) . 
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FIG. 1. The gluonic form factor F(p 2 ) as function of the momentum transfer (left panel: linear scale; right panel: log-log 
scale). Also shown is the solution of the set of DSEs proposed in [9] which have been solved for the case of SU(2) [10]. 

Physical units for the momentum can be obtained by using (|l5|). Calculations with different (3- values correspond to 
simulations with a different UV-cutoff A := 7r/a(/3). In order to obtain the renormalized gluon propagator, the gluonic 
wave function renormalization is chosen to yield a finite (given) value for the form factor at fixed momentum transfer. 

Figure [l] shows my result for the form factor F(p 2 ) where the condition ([ll] ) was implemented with an iteration over- 
relaxation method. The data obtained with different f3- values are identical within numerical accuracy, thus signaling 
a proper extrapolation to the continuum limit. 

At high momentum the lattice data are consistent with the behavior known from perturbation theory, 



iV) « i/ 



7 p2 

9 7 2 



13/22 



P 2 > 



(1 GeV) 2 



(19) 



The lattice data are compared with the solution of the gluon-ghost DSE |9|q From the DSE studies one expects a 
scaling law at small momentum 



F(p 2 ) oc [p 2 Y 



p 2 « (IGeV) 2 



(20) 



Depending on the truncation of the Dyson tower of equations and on the angular approximation of the momentum 
loop integral, one finds k = 0.92 || or re = 0.77 |ll| or re — > 1 Jl2]] . The lattice data are consistent with re = 0.5 
corresponding to an infrared screening by a gluonic mass. Also shown is the coarse grained "mass fit" (jj, = 1 GeV) 



Hp 2 ) 



p 



fi 



m\ [p 4 + m\ 



log 



13/22 



(21) 



* I thank Chr. Fischer for communicating his DSE solution for the SU(2) case prior to publication. 
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which nicely reproduces the lattice data within the statistical error bars. 



Results II: gluon propagator and confinement In order to get a handle on the information of quark confinement 
encoded in the gluon propagator in Landau gauge, I now change by hand the SU(2) Yang-Mills theory to a theory 
which does not confine quarks. It is instructive to compare the gluon propagator of the modified theory with the 
SU(2) result (see figure |). 

In the Maximal Center gauge ijjUf , the mechanism of quark confinement can be understood by a percolation of 
vortices which acquire physical relevance in the continuum limit Jl7[ . An intuitive picture in terms of vortex physics 
is also available for the deconfincmcnt phase transition at finite temperatures p8| . Reducing the full Yang-Mills 
configurations to their vortex content still yields the full string tension . Vice versa, removing these vortices from 
the Yang-Mills ensemble results in a vanishing string tension. This observation was used in []l9f to verify the impact 
of the vortices on chiral symmetry breaking. 
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FIG. 2. The static quark anti-quark potential (left panel) and the corresponding gluonic form factors (right panel); DSE 
solution from [10]. 



The static quark anti-quark potential in figure |^ demonstrates that a removal of the center vortices produces a non- 
confining theory. Figure [2] also shows the gluonic form factor obtained from the modified ensemble. The striking 
feature is that the strength of the form factor in the infra-red momentum range is drastically reduced. 

Results III: the Gribov noise Finally, let us check how strongly the gluonic form factor F(p 2 ) depends on the 
choice of gauge, i.e. on the sample of maxima of the variational condition (11) selected by the algorithm. For this 
purpose, I adopt an extreme point of view by comparing the gauge implemented by the iteration over-relaxation (IA) 
algorithm with the gauge obtained by simulated annealing (SA). The results of the form factor in both cases are 
shown in figure |[ I find, in agreement with [^3), that, in the case of the gluonic form factor, the Gribov noise is 
comparable with statistical noise for data generated with (3 6 [2.1,2.5] (scaling window). 
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FIG. 3. The gluonic form factor F(p 2 ) for a 10 3 x 20 lattice in the gauge IA and SA, respectively (left panel) and compared 
with previous results (IA, 16 3 x 32) (right panel). 
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